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^In  the  AURORA  flash  x-ray  machine,  a  bremsstrahlung  spectrum  is 
produced  in  four  thick  tantalum  targets  by  four  synchronous  10-MeV 
electron  beams.  The  bremsstrahlung  then  produces  an  electromag¬ 
netic  pulse  (EMP)  signal  through  the  same  mechanism  as  does  a 
nuclear  bomb.  It  induces  Compton  electron  currents  in  the  AURORA 
test  cell  chamber  (20  *  15  x  5  m).  As  the  Compton  electrons  are 
slowed  down,  they  ionize  the  air  and  also  produce  an  electromag¬ 
netic  field  in  the  chamber.  However,  the  Compton  current  is 
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20.  Abstract  (Cont'd) 

confined  to  a  relatively  small  volume  so  that  a  relatively  small 
electromagnetic  field  is  produced,  rendering  the  unmodified  AURORA 
machine  unsuitable  as  a  source-region  EMP  (SREMP)  simulator. 

A'  large  oarallel-plate  transmission  line  has  been  designed  to 
produce  a  r,rthreat-like"  EMP  pulse.  This  transmission-line  simu¬ 
lator  can  be  placed  in  the  AURORA  test  cell.  When  AURORA  and  the 
line  are  fired  synchronously,  an  SREMP-like  environment  is 
produced,  with  AURORA  providing  the  time-varying  air  conductivity. ! 

v-" 

The  configuration  described  above  has  been  used  to  study  the 
behavior  of  antennas  in  a  time-varying  conductive  medium.  Experi¬ 
mental  data  are  presented.  Also  presented  are  equivalent  circuit 
models  used  in  the  analysis  of  the  data.  Circuit  models  for 
antennas  are  obtained  using  the  "method  of  moments” — a  frequency- 
domain  technique.^  The  circuit  models  describe  the  antenna's 
interaction  with  the  environment  (through  the  equivalent  height 
circuit)  and  with  phe  loading  circuit  (through  the  impedance 
circuit).  The  modification  of  these  circuits  to  deal  with  time- 
varying  air  conductivity  is  discussed  (this  complicates  a 
frequency-domain  approach).  Comparisons  are  made  between  data 
taken  in  the  simulator  and  calculations  made  using  the  circuit 
models.  \ 
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INTRODUCTION 


Historically,  low-altitude  electromagnetic  pulse  (LAEMP)  research 
has  emphasized  strategic  interests.  Both  calculational  techniques 
(considering  only  a  localized  environment)  and  testing  (using  environ¬ 
ments  generated  during  underground  nuclear  tests)  have  been  widely 
applied  to  systems  expected  to  survive  extreme  overpressures.  Recently, 
the  possibility  o£  limited  nuclear  war  on  the  tactical  battlefield  has 
been  seriously  discussed,  and  interest  has  heightened  in  estimating  the 
response  of  tactical  equipment- -that  which  can  be  expected  to  survive  a 
much  less  severe  environment.  Unfortunately,  this  region  cannot  be 
effectively  simulated  with  underground  detonations.  Therefore,  atten¬ 
tion  has  been  directed  toward  other  methods  of  simulating  this  environ¬ 
ment  in  a  laboratory  setting,  as  well  as  developing  and  verifying  appli¬ 
cable  calculational  techniques. 

This  paper  describes  a  technique  for  the  simulation  of  the  tactical 
source-region  EMP  (SREMP)  environment  and  for  antenna  coupling 
analysis.  This  technique  is  being  developed  at  the  U.S.  Army's  Harry 
Diamond  Laboratories  (HDL) .  In  many  cases,  only  a  cursory  discussion  is 
provided,  with  many  of  the  technical  details  omitted.  Those  details 
will  be  published  in  a  series  of  technical  reports,  many  of  which  are 
presently  being  prepared. 


2.  TACTICAL  ENVIRONMENT 

It  has  been  recognized  that  no  tactical  system  can  be  expected  to 
survive  a  direct  nuclear  attack,  but  tactical  equipment  can  be  expected 
to  operate  near  nuclear  detonations.  A  method  has  been  developed  by  the 
U.S.  Army's  Nuclear  and  Chemical  Agency,  for  setting  attainable  surviv¬ 
ability  criteria  for  such  equipment.1  The  criteria  are  based  on  a 
"balanced  hardening"  philosophy  which,  simply  stated,  is  that  a  system 
should  be  hardened  against  any  particular  nuclear  effect  only  to  a  level 
commensurate  with  the  system's  hardness  against  all  other  effects.  For 
tactical  systems,  this  level  is  often  dominated  by  the  survivability  of 
the  system's  operator  to  blast,  thermal  radiation,  and  initial  radiation 
effects.  Systems  should  be  designed  to  survive  the  maximum  EMP  which 
could  be  experienced  at  a  range  (varying  for  different  yield  weapons) 
consistent  with  operator  survivability.  The  environment  is  character¬ 
ized  not  only  by  electromagnetic  fields  but  also  by  the  ionizing  radi¬ 
ation.  This  ionizing  radiation  complicates  the  electromagnetic  effects 
by  introducing  current  sources  and  time-varying  air  conductivity. 


Nuclear  Survivability  Criteria  for  Army  Tactical  Equipment  (U),  U.S . 
Army  Nuclear  Agency,  NUA-C-5400-74  (August  1974).  (CONFIDENTIAL- 
RESTRICTED  DATA) 


3.  AURORA  FLASH  X-RAY  MACHINE 


The  Defense  Nuclear  Agency's  AURORA  Flash  X-Ray  Facility  at  HDL  in 
Adelphi,  MD,2  **  is  the  largest  of  its  type.  This  machine  produces 

bremsstrahlung  in  four  thick  tantalum  targets  by  four  synchronous  pulsed 
8-MeV,  230-kA  electron  beams.  At  a  total  energy  of  about  1.2  MJ,  the 
machine  was  designed  to  provide  a  high-dose  gamma  radiation  environment 
over  approximately  a  1-m  cube  (shown  later,  sect.  4.1,  fig.  4(a)). 

Fortunately,  the  facility  is  equipped  with  a  large  test  cell  (19.8  x 

14.6  x  4,6  m),  rendering  it  suitable  for  experiments  requiring  lower 
levels  of  radiation  over  a  considerably  larger  volume.  The  radiation 
and  resulting  Compton  currents,  acting  through  the  same  physical 

mechanisms  as  expected  from  a  nuclear  detonation,  produce  an  EMP 

environment  in  the  facility's  test  chamber  which  is  qualitatively 
similar  to  that  resulting  from  a  nuclear  detonation.  However,  because 
of  several  factors,  including  the  time  history  of  the  radiation,  the 
geometry  of  the  sources,  the  small  volume  of  Compton  current,  and  the 
conductive  walls  of  the  test  chamber,  the  unmodified  AURORA  facility  is 
unsuitable  as  a  source-region  simulator.  Figure  1  shows  a  typical 

normalized  gamma  radiation  pulse  generated  by  AURORA.  The  resulting  air 
conductivity  waveform  has  approximately  the  same  shape,  with  peaks 

varying  from  10”1 *  mho/m  at  the  "hot  spot"  to  about  10-1*  mho/m  at  the 

back  of  the  test  chamber.  This  gamma  pulse  time  history  differs  from 
the  prompt  radiation  which  would  produce  the  survivability  criteria 
environment. 

At  a  point  in  the  chamber  where  the  peak  conductivity  is  approxi¬ 
mately  that  associated  with  operator  survivability  (10  m  from  the 
source),  the  time  histories  of  the  conductivity  and  vertical  electric 
fields  differ  significantly  from  those  of  the  survivability  criteria. 
Figure  2  shows  the  vertical  electric  field  and  conductivity  at  the  point 
as  predicted  by  AUR3D,  a  computer  code  produced  by  William  Crevier  of 
Mission  Research  Corporation  to  simulate  the  AURORA  environment.  These 
results  compare  favorably  with  experimental  results  obtained  by  HDL  at 
the  same  location.5  Neither  the  electric  field  nor  the  conductivity 
produced  by  AURORA  model  well  the  electric  field  and  conductivity  which 
are  expected  in  a  typical  tactical  EMP  environment.  They  differ  sub¬ 
stantially  in  peak  amplitude,  risetime,  and  late-time  behavior. 


2Hans  Fleischmann ,  High  Current  Electron  Beams,  Physics  Today,  28,  5 
(May  1975),  34. 

*Gerold  Yonas,  Fusion  Power  with  Particle  Beams ,  Scientific  American, 
(November  1978),  50. 

**B.  Bernstein  and  I.  Smith,  AURORA,  An  Electron  Accelerator,  IEEE 
Trans  Nucl  Sci,  NS-20 ,  3  (June  1973). 

5L.  M.  Anderson  et  al ^  Status  of  AURORA  Environment  Calculations, 
IEEE  Trans  Nucl  Sci,  NS-2 5 ,  6  (December  1978). 


TIME  AFTER  ZERO  TIME 


Figure  1*  AURORA  gamma  pulse.  A  typical  pulse 
can  be  expected  to  fall  within  the  shaded  area 


Figure  2.  E-field,  its  time  derivative,  and  air 
conductivity  at  base  of  antenna  at  10  m  as  calculated 
by  AUR3D  with  no  antenna  present. 


Consideration  of  a  short  monopole  antenna  demonstrates  the  inade¬ 
quacy  of  the  AURORA  environment. 6  The  short-circuit  current  in  a  short 
antenna  is  given  as 

_  1  ../dE  ,  oE\ 

xs  =  7  hc  \Xt  +  T^)  ’ 

where 

h  =  the  height  of  the  antenna 
C  =  the  capacitance  of  the  antenna  =  —s— 

ZQ  =  60[2Jln[|^)  -  3.39“  ° 

a  =  the  antenna  wire  radius 
E  =  the  incident  electric  field 
o  =  the  conductivity  of  the  medium 
e0  =  the  permittivity  of  the  medium 
c  =  velocity  of  light 

Figure  3  shows  the  response  of  the  short  monopole  antenna  exposed  to 
the  AURORA  environment  at  a  point  where  the  peak  conductivity  is  consis¬ 
tent  with  survivability  criteria.  The  good  agreement  between  AURORA 
data  and  this  model,  which  does  not  include  any  resonant  behavior, 
suggests  that  the  unmodified  AURORA  is  unsuitable  as  an  SREMP 
simulator.  In  a  typical  tactical  environment,  resonances  would  be 
excited  on  a  1-m  monopole  antenna  such  as  the  one  used  in  gathering 
these  data. 

4.  MODIFICATION  OF  AURORA  ENVIRONMENT 
4.1  First-Order  Modification 

It  is  apparent  from  the  foregoing  discussion  that  E-fields, 
which  are  larger  and  have  shorter  risetimes  than  those  available  from 
unmodified  AURORA,  are  essential  for  even  the  most  rudimentary  simul¬ 
ation.  However,  AURORA  is  the  only  available  facility  which  can  provide 
acceptable  time-varying  conductivity  levels  over  a  large  enough  volume 
to  be  useful.  It  appears  that,  in  order  to  take  advantage  of  this 
unique  feature  of  AURORA,  one  should  make  use  of  a  supplementary  EMP 
source.  For  the  work  reported  here,  this  supplementary  source  took  the 
form  of  a  large  (12  x  3  x  5  m)  parallel -plate  transmission  line — 
designed  with  the  assistance  of  William  Crevier.  The  two  plates  of  this 
line  were  mounted  vertically  in  the  test  chamber  with  a  separation  of 

6fi.  Manriquez ,  G.  Merkel,  W,  D.  Scharf ,  and  D.  Spohn,  Electrically 
Short  Monopole  Antenna  Response  in  an  Ionized  Air  Environment— 
Determination  of  Ionized  Air  Conductivity ,  IEEE  Trans  Nucl  Sci,  26 ,  6 
(December  1979), 


11 


•  •  • 


3  in  and  stretched  laterally  across  the  chamber  from  one  side  wall  to  the 
other.  (See  fig.  4(a)  to  (c).)  The  hindmost  plate  (furthest  from  the 
"hot  spot")  is  grounded  to  the  test  cell  floor  and  ceiling,  but  the 
front  plate  is  isolated  from  ground  and  is  driven  by  a  100-kV  pulser. 
The  distance  between  the  hot  plate  and  both  the  floor  and  ceiling  is 
about  1  m,  and  the  plate  curves  in  at  both  the  top  and  bottom,  so  that 
the  plate  separation  at  top  and  bottom  is  only  about  230  cm.  This  cur¬ 
vature  focuses  the  field  lines  so  that  a  relatively  uniform  E-field  is 
achieved  at  half-height,  where  experimental  apparatus  is  mounted.  In  an 
attempt  to  minimize  impedance  mismatches,  both  the  ground  plate  and  the 
hot  plate  were  tapered  out  from  the  drive  point.  The  line  tapers  out 
from  a  separation  of  10  cm  at  the  drive  point  to  the  full  3-m  separation 
over  a  distance  of  4  m.  The  line  was  terminated  at  the  other  end  in  an 
anechoic  resistive  curtain  connecting  the  two  plates  and  in  two  resis¬ 
tive  bundles  connecting  the  hot  plate  to  the  floor  and  ceiling.  The 
termination  was  designed  to  distribute  the  current  roughly  along  field 
lines.  The  entire  transmission  line  was  constructed  as  a  unit  on  a 
sturdy  wooden  frame,  so  that  it  could  be  moved  back  and  forth  in  the 
test  cell — providing  a  means  of  adjusting  the  conductivity  level.  The 
line  was  driven  at  the  drive  point,  through  50  ft  of  50-ohm  cable,  from 
a  100-kV  pulser  built  by  Pulsar,  Inc.,  which  delivers  a  pulse  with  a 
risetime  of  about  20  ns  and  a  decay  time  of  about  2  us.  Thus,  a  peak 
field  of  about  33  kV/m  is  seen  in  the  line. 

It  should  be  noted  that  the  E-field  signal  is  by  no  means 
independent  of  the  AURORA-induced  air  conductivity.  In  order  to  analyze 
the  effect  of  ionizing  radiation  on  the  transmission  line,  we  must  know 
the  air  conductivity  as  a  function  of  ionizing  radiation.  In  1977,  HDL 
carried  out  a  number  of  "pie-pan"  experiments7  to  measure  air  conduc¬ 
tivities  for  a  range  of  humidities.  This  work,  analyzed  by  Dean  May,8 
is  summarized  in  figure  5. 

As  an  aid  to  visualizing  the  effect  of  the  conductivity, 
consider  the  simple  circuit  model  of  figure  6.  (This  is  intended  only 
as  a  very  rough  qualitative  description  of  the  interaction.)  In  this 
model,  the  pulse  (which  has  a  rlong  discharge  time  «  2  us)  is  assumed  to 
be  a  perfect  voltage  source.  The  50-ohm  cable  is  modelled  as  an 
inductor,  the  line  as  a  capacitor,  and  the  termination  as  a  resistor. 
(Of  course,  in  reality,  both  the  line  and  the  cable  have  distributed 
inductance  and  distributed  capacitance.) 


7W.  F,  Crevier,  C.  L.  Longmire,  G.  Merkel,  D.  J.  Spohn,  Air  Chemistry 
and  Boundary  Layer  Studies  with  AURORA,  IEEE  Trans  Nucl  Sci,  MS-2 4 ,  6 
(December  1977), 

8Dean  L,  May,  Electron  Mobility  and  Electron  Attachment  Rate  in  an 
Electromagnetic  Pulse  Air „  Environment ,  Harry  Diamond  Laboratories 
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Figure  5.  Theoretical  and  experimental  ratios  of  electron 
mobility  to  electron  attachment  rate  versus  E-field  for  dry  and 
moist  air.  To  a  very  good  approximation,  air  conductivity  is 
given  by  o(E)  -  yqu(E)/a(E),  where  y  is  the  ionization  rate  in 
ion  pairs  per  second.  The  3.8-percent  water  vapor  curve  is 
nearly  flat;  i.e.,  conductivity  at  3.8-percent  water  vapor  is 
independent  of  electric  field.  Therefore,  to  a  very  good 
approximation,  at  3.8-percent  water  vapor,  air  conductivity, 
o(t),  is  given  by  o(t)  =  0.5  x  10“12  Fy(t),  where  Fy  is  in 
roentgens  per  second.  Theoretical  curves  are  calculated 
in  the  manner  described  in  Air  Chemistry  and  Boundary  Layer 
Studies  with  AURORA,  by  Crevier,  Longmire,  Merkel,  and  Spohn, 
IEEE  Trans.  Nucl.  Sci.,  NS-24  (December  1977). 
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Figure  6.  Simplified  equivalent  circuit  of 
transmission- line  simulator. 


If  the  capacitor  dielectric  (air)  is  subjected  to  pulsed 
ionizing  radiation  causing  a  peak  conductivity  of  10-3  mho/m,  the  capa¬ 
citor  "turns  into"  a  6-ohm  resistor.  When  the  conductivity  disappears 
(when  the  radiation  stops)/  it  turns  back  into  a  capacitor. 
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The  circuit  appearing  in  figure  7  can  therefore  be  described 
even  more  simply  by  the  following  equation: 


V  =  L  — ~  +  RI  . 
at 


This  equation,  in  which  V  and  R  are  functions  of  time,  can  be 
solved  using  the  method  of  variation  of  parameters: 


I(t)  =  e 


-IS  *<*>«*  . JS  ‘ft*" 


In  v<*> 


dt 


where 


R(t)  .  /fc.  V(t) 
r(t)  -  — and  v(t)  =  — ~ 


Pigure  7.  Extremely  simplified 
equivalent  circuit  of  transmission¬ 
line  simulator. 


Qualitatively,  the  reaction  of  the  circuit  to  a  pulsed  dip  in 
load  resistance,  with  V  constant,  can  be  described  as  follows.  As  the 
resistance  falls,  a  magnetic  field  (and  hence  a  voltage  drop)  builds  up 
in  the  inductor.  The  result  is  a  drop  in  the  voltage  across  the  load. 
After  the  current  reaches  its  peak,  the  energy  stored  in  the  magnetic 
field  discharges  back  into  the  circuit,  the  current  rises,  and,  if  the 
load  resistance  has  risen  back  to  its  full  value  in  a  time  small 
compared  to  the  discharge  time  of  the  inductor  (L/R,  where  R  is  an 
effective  "average"  value  of  the  resistance),  the  voltage  across  the 
load  will  overshoot  its  quiescent  level.  This  "inductive  kick"^'10  can 
be  seen  in  both  the  measured  and  the  calculated  transmission-line 
voltage . 

4.2  Suggested  Further  Modifications 

The  above  description  applies  to  the  case  of  a  constant  voltage 
source.  If  there  is  appreciable  conductivity  during  the  risetime  of  the 
pulse,  the  desired  pulse  characteristics  can  be  seriously  degraded.  At 
the  criterion  level  of  3  x  i  o*4  mho/m  peak  conductivity,  this  degrada¬ 
tion  begins  to  be  a  problem.  Figure  8  compares  the  E-field  measured  in 
the  transmission  line  when  the  AURORA  is  fired  to  the  E-field  when  the 
AURORA  is  not  fired.  Figure  9  is  a  similar  comparison,  but  with  the 
line  3  m  closer  to  the  AURORA  hot  spot.  At  this  spot,  the  conductivity 
is  much  higher  and  there  is  a  correspondingly  greater  E-field  distor¬ 
tion.  (Note  that  the  relative  firing  times  of  AURORA  and  the  pulser 
differ  in  the  two  figures.  AURORA  fired  first  in  fig.  8,  while  the 
pulser  fired  first  in  fig.  9.)  Since  measurements  are  desired  in  an 
environment  characterized  both  by  conductivity  and  fast-rising  fields, 
steps  must  be  taken  to  minimize  this  waveform  degradation.  There  are  a 
number  of  approaches  to  the  problem  of  reducing  the  inductive  kick.  The 
most  obvious  method  is  to  position  the  line  so  that  the  conductivity  is 
low.  This  clearly  limits  the  utility  of  the  transmission-line  approach, 
since  time-varying  conductivity  is  an  integral  part  of  the  environment 
to  be  simulated. 

The  inductive  kick  can  also  be  limited  by  reducing  the  source 
impedance  of  the  pulser  (e.g.,  raising  the  capacitance)  and  shortening 
the  50-ohm  feeder  cable  (lowering  its  inductance).  However,  one  is 
limited  ultimately  by  the  distributed  inductance  of  the  transmission 
line  itself.  This  can  be  handled  by  introducing  a  "slave  line,"  which 
shares  the  same  hot  plate,  but  which  is  filled  with  a  nonionizing  gas 


Manriquez ,  G.  Merkel,  and  D.  Spohn,  Modification  of  the  AURORA 
Environment,  Harry  Diamond  Laboratories,  HDL-PR-79-5  (October  1979). 

1  °R.  Manriquez ,  G .  Merkel,  and  D .  Spohn,  Modification  of  the  AURORA 
Electromagnetic  Environment ,  IEEE  Annual  Conference  on  Nuclear  and  Space 
Radiation  Effects  (July  1978). 


such  as  SFfc  or  CO2.  Physically,  this  can  be  achieved  by  erecting 
another  ground  plate  in  front  of  the  hot  plate  and  filling  it  with  the 
appropriate  gas  (see  fig.  10).  The  resulting  line  can  then  be  visual¬ 
ized  as  a  "coaxial  cable"  filled  with  nonionizing  gas,  except  for  an 
"experimentation  sector,"  which  is  filled  with  ionizing  gas  (such  as  air 
or  even  N2  for  higher  conductivity.  See  fig  11.)  Of  course,  this 
"master-slave  line"  (or  "sectored  coax")  will  have  a  much  lower 
impedance  than  the  simple  vertical-plate  line.  Thus,  it  will  "waste" 
energy  and  make  more  demands  on  the  pulser.  The  pulser  impedance 
becomes  a  critical  consideration.  HDL  has  a  variable  voltage  (40  to  400 
kV)  pulser  with  a  decay  time  of  1.5  11s  into  10  ohms.  This  pulser  should 
be  suitable  for  use  with  the  sectored  line. 

It  should  also  be  mentioned  that  the  ionization  enhancement 
techniques  outlined  above  may  have  implications  for  the  study  and  simu¬ 
lation  of  the  SREMP  environment  in  the  strategic  case — for  example,  the 
MX  System  or  the  Ballistic  Missile  Defense  System.  Ordinarily,  such 
systems  (or  their  subsystems)  could  be  tested  only  in  an  underground 
test.  With  our  proposed  methods,  a  much  more  cost-effective  and  repeat- 
able  series  of  tests  could  be  conducted,  which  would  reliably  simulate 
at  least  the  purely  electromagnetic  effects  on  such  strategic  systems. 
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Figure  8.  Comparison  of  E-field  with  and  without  ambient  time-varying 
ionized  air  conductivity. 
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Figure  9.  Comparison  of  E-field  with  and  without  ambient  time-varying 
ionized  air  conductivity  (3  m  closer  to  hot  spot  than  in  figure  8.) 


In  the  application  of  the  scheme  described  above  to  strategic 
source-region  simulation,  there  will  be,  for  very  high  air  conductivity, 
a  self-shielding  effect  which  will  inhibit  the  propagation  of  the  trans¬ 
mission-line  mode.  This  effect  becomes  critical  when  the  diffusion 
length  over  the  time  of  interst  is  of  the  same  order  of  magnitude  as  the 
transmission-line  dimensions.  Then,  if  the  time  of  interest  is  100  ns 
and  a  typical  distance  is  5  m. 
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Figure  10.  Comparison  of  present  and  proposed  new  improved 
transmission  lines  that  would  not  suffer  the  "inductive  kick' 
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Figure  11.  Generalization  and  simplification  of  slave  line 
concept.  The  smaller  the  "work"  section,  the  less  the 
"inductive  kick." 

Thus,  as  long  as  the  conductivity  remains  below  about  0.01  mho/m,  the 
transmission  line  should  operate  as  described.  For  higher  conductivity, 
the  fields  are  determined  (as  they  are  in  true  nuclear  EMP)  by  the  local 
Compton  current — i.e.,  the  current  within  a  diffusion  sphere.  As  time 
progresses,  this  sphere  ( j R  -  RQ  |  <  lyiff)  expands  for  two  reasons: 

(1)  the  explicit  square-root  dependence  of  on  t  and  (2)  the 

implicit  dependence  through  o(t),  which  decreases  at  late  time. 

For  high  conductivity  levels,  where  the  line  self-shields,  the 
most  promising  method  of  simulation  seems  to  be  the  direct  injection  of 
high-energy  electrons  into  the  test  volume.  For  instance,  AURORA  could 
be  operated  in  the  electron-beam  mode,  which  is,  in  any  case,  a  far  more 
efficient  use  of  the  pulsed  energy  than  is  achieved  in  the  more  conven¬ 
tional  bremsstrahlung  mode.  Electron  irradiation  of  a  large  volume  can 
be  achieved  through  the  use  of  a  high-Z  thin-target  multiple-scattering 
"spray  nozzle." 

It  should  be  noted  that  the  proposed  use  of  nonattaching  gases 
will  increase  the  fidelity  of  the  o(t)  waveform  in  two  ways:  (1)  by 
causing  a  levelling-off  arising  from  the  constant  resupply,  through 
Compton  creation  and  slowing  down,  of  conduction  electrons  and  (2)  by 
enhancing  the  rise  rate  of  the  conductivity  waveform.  The  free  electron 
density  n(t)  obeys  the  equation 


a* 


Q  +  an 


where  Q  is  the  ion-pair  creation  rate,  and  a  is  the  attachment  rate. 
(More  generally,  other  recombination  terms  may  be  included.)  The  exact 
solution,  with  initial  condition  n(0)  *  0,  is 

n(t)  =  exp  [-at]  exp[at] Q( T)dt  . 


Note  that,  although  the  risetime  is  longer  for  low  attachment  rates,  the 
rise  rate  is  greater,  leading  to  a  greater  contribution  from  the 
parametric-drive  e)  term.  For  xenon,  the  a  =  0  curve  applies,  and 
thus  has  the  same  shape  as  that  for  nitrogen,  but  the  magnitude  is 
enhanced  by  a  factor  of  5.5,  due  to  higher  Z  (electron  density)  and  due 
to  a  smaller  electron  energy  loss  per  ion  pair  (36.3  eV  for  N2,  21.4  eV 
for  Xe) . 


In  choosing  gases  for  use  in  SREMP  simulation,  one  must 
consider  carefully  their  behavior  when  ionized.  One  important  effect  is 
the  formation  of  boundary  layers.  In  a  low-density  plasma,  the  well- 
known  Debye  sheath  effect  occurs  at  a  metallic  surface.  In  a  plasma  at 
atmospheric  pressure  (collision-dominated),  the  layer  formation  is 
dominated  by  a  number  of  medium-dependent  parameters — e.g.,  plasma 
frequency,  collision  frequency,  ionic  and  electronic  recombination  and 
attachment  rates,11  ionic  and  electronic  mobilities,12  avalanching 
properties,  etc.  Electrode  geometry  is  also  very  important.  Small 
radii  of  Curvature  intensify  fields  and  thus  intensify  both  boundary- 
layer  formation  (through  electron  depletion)  and  avalanching.  These  two 
effects  oppose  each  other  so  that,  for  instance,  in  nitrogen  gas,  the 
use  of  grid  electrodes  actually  reduces  the  boundary-layer  effects,  as 
compared  to  the  use  of  plate  electrodes.  These  phenomena,  whose 
descriptions  here  have  been  a  bit  oversimplified,  are  currently  under 
study  at  HDL,  and  will  be  reported  on  soon.*  Figure  12(a)  compares  the 
current  through  ionized  air  and  through  nitrogen  around  grid  and  plate 
electrodes  in  a  large  biased  capacitor  exposed  to  ionizing  radiation. 
Figure  12(b)  shows  the  apparatus  used  to  obtain  the  data  given  in  figure 
12(a). 


1 1H.  J.  Longle y  and  C.  L.  Longmire,  Electron  Mobility  and  Attachment 
Rate  in  Moist  Air,  Mission  Research  Corporation,  MRC-N-222  (1975). 

12C.  E .  Baum,  The  Calculation  of  Conduction  Electron  Parameters  in 
Ionized  Air,  Electromagnetic  Pulse  Theoretical  Notes,  I_,  Note  6  (May 
1967). 

*M.  Bushel  1 ,  C .  Kenyon,  G.  Merkel,  and  W.  D .  Scharf ,  Electron 
Depletion  Effects  in  the  Nuclear  EMP  Source  Region.  To  be  submitted  to 
IEEE  Trans  Nucl  Sci. 
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Figure  12.  Influence  of  grid 
structure  on  "pie-pan"  conductivity 
measurements:  (a)  Peak  current  as 
a  function  of  bias  voltage  in  air 
and  nitrogen. 
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Figure  12  (cont'd).  Influence  of  grid  structure  on  "pie-pan" 
conductivity  measurements:  (b)  Experimental  setup  used  to  obtain  data 
shown  in  figure  12(a) «  Data  marked  "center  grid"  in  figure  12(a)  were  9 

obtained  with  grid  in  "pie-pan"  rather  than  plate.  Large  value  of  —  •  - 

nitrogen  data  with  grid  in  center  of  "pie-pan"  results  from  avalanching 
breakdown  within  the  electron  depletion  boundary  layer. 
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5.  THEORETICAL  INTERPRETATION  OF  INTERACTION  OF  TRANSMISSION-LINE  AND 
IONIZING  RADIATION 

The  detailed  behavior  of  the  transmission  line  when  the  AURORA  is 
fired  has  been  interpreted  with  a  36-section  lumped -para me ter  trans- 
mission-line  model.1 3  The  first  18  sections  have  been  used  to  model  the 
50-ft  50-ohm  coaxial  cable,  and  the  1 2-m-long  parallel  plate  has  been 
modelled  with  a  corresponding  18-section  lumped-parameter  network  (LPN) 
transmission  line.  The  best  preliminary  fit  to  the  measured  data  was 
obtained  when  the  parallel-plate  line  was  assumed  to  be  about  65  ohms. 
The  line  was  terminated  with  a  67-ohm  resistive  screen.  Figure  13  is  a 
schematic  diagram  of  the  parallel-plate  LPN  model;  figure  14  compares 
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Figure  13.  Schematic  of  transmission-line  lumped  parameter  network 
model. 


13*.  Bushell,  R.  Manriquez,  G.  Merkel,  W,  Scharf ,  and  D.  Spohn, 
Source-Region  EMP  Simulator — A  Parallel  Plate  Transmission  Line  in  the 
AURORA  Test  Cell,  IEEE  Trans  Nucl  Sci,  MS-27,  6  (December  1980). 


the  E-field  measured  inside  the  parallel-plate  transmission  line  and  the 
value  calculated  with  our  LPN  model  when  no  conductivity  is  present. 


A  time-domain  reflectometry  { TDR)  analysis  of  the  transmission  line 
indicates  that  the  transition  region  between  the  coaxial  cable  and  the 
parallel -plate  line  (the  tapered  section)  has  a  varying  impedance.  In  a 
forthcoming  analysis,  we  will  attempt  to  incorporate  this  information 
into  our  analysis,  and  we  hope  to  reproduce  the  "fine  structure"  of  the 
measured  data. 
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Figure  14.  E-field  inside  parallel-plate  transmission  line  (no 
conductivity) . 
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The  effect  of  ionizing  radiation  on  the  parallel -plate  line  is 
incorporated  into  the  line  model  by  introducing  a  time-varying  resistor 
across  each  capacitor.  Hie  relationship  between  the  value  of  this 
resistance  and  the  ambient  conductivity,  o(t),  is  given  by 

R(t)  =  — — —  . 

C  a(t) 

This  equation  is  a  direct  consequence  of  the  fact  that  the  value  of  the 
capacitance  and  the  value  of  the  resistance  obtained  by  filling  the 
capacitor  with  a  conductance  are  both  obtained  by  solving  Laplace's 
equation. 

Figures  15  to  20  compare  an  experimental  E-field  measurement  and 
results  calculated  from  the  LPN  model.  The  same  experimental  curve 
appears  on  all  six  figures.  The  theoretical  curve  of  figure  15  uses  a 
conductivity  time  history  obtained  with  radiation  dosimetry  measurements 
and  with  the  air  chemistry  parameters  given  by  Dean  May's  analysis. 
Figures  16  to  20  give  theoretical  results  using  renormalizations  of  the 
same  a(t)  curve,  each  lower  than  the  previous  one  by  a  factor  of  two. 
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Figure  17.  E-field  inside  parallel-plate  transmission  line  f a___v  =  1  * 
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Figure  20.  E-field  inside  parallel-plate  transmission  line  f o„alr  = 
1.25  x  10"4  mho/m).  136 


6.  ANTENNA  COUPLING  EXPERIMENTS  AND  ANALYSIS 


6.1  Experimental  Procedures 


In  a  typical  test  shot  procedure,  an  antenna  short-circuit 
current  response  was  first  measured  when  the  parallel -plate  line  was 
fired,  without  firing  the  AURORA.  Next,  the  measurement  was  repeated 
with  AURORA  and  the  line  firing  simultaneously.  As  already  indicated, 
figure  15  shows  a  typical  E-field  measurement  obtained  with  a  parallel- 
plate  sensor  whose  bottom  plate  is  the  transmission-line  ground  plate 
itself  and  whose  top  plate  is  a  round  sheet  of  aluminum  foil  positioned 
at  a  1/2-in.  separation  from  the  ground  plate.  This  sensor  has  a  capa¬ 
citance  of  80  pF  and  is  loaded  with  a  10-kohm  resistor.  E-field  sensors 
were  placed  at  three  positions  in  the  transmission  line.  Hie  magnetic 
field  in  the  parallel-plate  transmission  line  was  also  measured  at  the 
same  three  positions  with  Moebius  loop  magnetic-field  sensors.  The 
current  flow  into  the  transmission  line  was  also  monitored.  Unfortun¬ 
ately,  because  of  the  large  magnitude  of  current  involved,  the  current 
probe  did  not  perform  up  to  expectations.  A  photodiode  radiation 
detector  was  placed  under  the  transition  from  the  coaxial  cable  to  the 
parallel-plate  transmission  line  (see  fig.  4(a)).  The  photodiode 
yielded  an  absolute  value  of  the  radiation  dose  rate  in  rads/s.  Thirty 
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to  forty  thermoluminescent  radiation  detectors  (TLD's)14  were  placed 
throughout  the  test  cell  to  monitor  the  total  radiation  dose  at  selected 
points  throughout  the  room.  The  transmission  line  was  terminated  by  a 
matched  load  consisting  of  a  curtain  of  parallel  resistor-loaded 
wires.  The  current  through  the  curtain  was  also  measured. 
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To  obtain  the  antenna  coupling  data  presented  here,  linear 
monopole  antennas  of  various  lengths  were  mounted  horizontally  at  half¬ 
height  in  the  transmission  line  and  shorted  to  the  ground  plate.  The 
radiation  dose  was  monitored  by  four  TLD's  attached  to  the  antenna,  one 
at  the  tip,  one  at  the  base,  and  two  intermediate.  Current  response 
measurements  were  made  with  a  Singer  probe  (transfer  impedance  -  0.1 
ohm)  mounted  around  the  base  of  the  antenna.  The  effect  of  the  ionized 
air  on  the  interaction  of  the  antenna  with  an  incident  E-field  can  be 
summarized  by  saying  that  the  air  conductivity  tends  to  cause  bipolar 
(oscillatory)  responses  to  become  monopolar  as  a  result  of  the  dominance 
of  conductive  current  (oE)  over  displacement  current  (e  ) .  The  oscil¬ 
lations  also  tend  to  become  damped.  The  detailed  theoretical  expla¬ 
nation  will  be  presented  in  later  sections. 

6.2  Equivalent  Circuit  Method 

Hie  authors  have  developed  a  method  of  modelling  antenna 
coupling  in  time-varying  conductive  air.  This  method  uses  equivalent 
circuit  models  obtained  by  established  methods.15  Hiese  circuits  are 
modified  according  to  the  algorithm  described  in  this  section. 

Often,  the  analysis  of  antennas  is  undertaken  in  the  frequency 
domain.  This  is  a  natural  course  of  action  for  dealing  with  systems 
whose  properties  (e,  u,  o)  do  not  change  with  time,  for  in  such  a  case 
all  differential  (or  integro-differential)  equations  which  arise  have 
time-independent  coefficients  and  kernels.  Fourier-transforming  these 
equations  converts  all  time-differential  or  time-integral  operators  into 
simple  multipliers,  greatly  facilitating  the  analysis.  If,  in  addition 
to  being  time-independent,  the  material  properties  are  also  space- 
independent,  the  difficulty  of  the  analysis  is  determined  only  by  the 
system  geometry. 

A  frequency-domain  analysis  typically  seeks  to  describe  an 
antenna  by  two  complex  functions  of  frequency — h(a>)  and  Z(oi).  The 
former  is  the  effective  height,  and  describes  the  antenna's  "cross 
section"  (its  interaction  with  its  environment).  (This  environment 
includes  the  spatial  variation  of  the  electromagnetic  driving  signal  so 


^Klaus  G •  Kerris,  The  AURORA  Dosimetry  System,  Harry  Diamond 
Laboratories,  HDL-TR-1754  (March  1976). 

15John  Sweton,  DMA  EMP  Handbook,  Vol .  2,  Coupling  Analysis  (U), 
Defense  Nuclear  Agency ,  DNA  2114H-2 ,  Ch  10  (July  1979). 
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that  the  effective  height  would  depend,  for  example,  on  the  angle  of 
incidence  of  a  plane-wave  signal.)  The  latter  is  the  antenna  impedance, 
and  describes  the  antenna's  behavior  as  a  circuit  element — its 
interaction  with  its  loading  circuit. 

A  variety  of  techniques  exist  for  obtaining  the  effective 
height  and  impedance  functions.  Among  these  is  the  "method  of  the 
finite-element"16'17  used  by  mechanical  engineers.  The  authors  have 
chosen  to  use  this  method  in  their  analysis.  The  effective  height  and 
impedance  circuits  for  a  vertical  monopole — valid  for  frequencies  around 
and  below  the  first  resonance--are  shown  in  figures  21  and  22.  These 
circuit  models  are  consistent  with  those  presented  by  P.  P.  Toulios. 18 

CS 


Figure  21.  Antenna  effective  length  equivalent  circuit  for 
first  antenna  resonance.  Effect  of  time-varying  ambient 
conductivity  on  antenna  effective  height  is  incorporated 
into  above  circuit  by  inserting  time-varying  resistor  in 
parallel  with  each  of  the  three  capacitors  CA(i  =  1,2,3): 
RA(t)  =  e/[o(t)Cil.  Equivalent  circuits  of  this  type  are 
not  limited  to  simple  linear  antennas. 


1 ^R.  F\  Harrington,  Field  Computation  by  Moment  Methods,  New  York, 
Macmillan  (1968). 

17Hu.  H.  Chao  and  Bradley  J.  Strait,  Computer  Programs  for  Radiation 
and  Scattering  by  Arbitrary  Configurations  of  Bent  Wires,  Electrical 
Engineering  Dept.,  Syracuse  University ,  Air  Force  Cambridge  Research 
Laboratories  AFCRL-7 0-0  37  4  ( September  1970).  (The  Fortran  code 

described  by  Chao  and  Strait  has  been  altered  so  that  it  can  be  used  for 
a  medium  with  constant  conductivity . ) 

18P.  P.  Toulios,  Antenna  User's  Manual  for  Linear  Cylindrical  Antennas 
in  an  EMP  Environment,  Vol .  I,  IIT  Research  Institute,  DAAG39-72-C-0 192 
(January  1974). 
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Figure  22.  Antenna  impedance  equivalent  circuit  for  first 
antenna  resonance.  Effect  of  time-varying  ambient 
conductivity  on  antenna  impedance  is  incorporated  into 
above  circuit  by  inserting  time-varying  resistor  in 
parallel  with  each  of  the  two  capacitors  C^(i  =  1,2):  Ri  ( t ) 
=  e/[o(t)CjJ.  Equivalent  circuits  of  this  type  are  not 
limited  to  simple  linear  antennas. 


Once  h(w)  and  Z(w)  are  found,  it  is  useful  to  synthesize 
passive  circuits  which  model  the  antenna.  These  circuits  serve  two 
purposes:  (1)  to  aid  physical  intuition  and  (2)  to  enable  workers  to 

use  standard  network  analysis  codes  in  predicting  antenna  response.  An 
antenna  equivalent  circuit  can  easily  be  hooked  onto,  say,  a  radio 
circuit  to  determine  the  total  system  response.  Equivalent  circuit 
models  have  been  obtained  by  Toulios  for  a  number  of  antennas  used  with 
Army  equipment.  These  models  have  proven  very  useful  in  dealing  with 
high-altitude  EMP  signals,  but  are  clearly  inappropriate  for  the  source- 
region  case,  for  which  time-varying  air  conductivity  must  be  taken  into 
account. 


Any  antenna  has  a  dc  capacitance  and  an  infinite  set  of 
resonances.  It  seems  intuitively  natural  to  represent  this  as  a  Foster 
canonical  circuit. 19,20  Of  course,  such  a  representation  can  only 
describe  a  purely  reactive  system,  and  must  be  modified  (by  the  addition 
of  resistors)  to  include  the  effects  of  radiation  resistance  and  of 
finite  air  conductivity.  (See  fig.  23.  We  ignore  antenna  and  ground 

191*.  E.  Van  Valkenburg ,  Introduction  to  Modern  Network  Synthesis ,  Ch. 
5,  New  York,  John  Wiley  and  Sons,  Inc,  (1967), 

20R.  M,  Foster,  A  Reactance  Theorem,  Bell  System  Technical  Journal,  3_ 
(April  1924),  259-267. 


Figure  23.  First  Foster  form  of  reactive  network.  CQ  is  present  if 
Z(w)  has  pole  at  o>  **  0,  and  L—  is  present  if  Z (o>)  has  pole  at  infinity. 
After  Van  valkenburg  (ref.  19). 


resistivity.)  Ihe  resistors  which  represent  air  conductivity  appear  in 
parallel  with  each  capacitor  in  the  model  in  accordance  with  the  pre¬ 
scriptive  formula:  R  =*  e0/oC.  Thus,  a  purely  capacitive  admittance  ( Y 
*  ioC)  acquires  a  dissipative  term  (y  *  iwC  +  -^).  This  corresponds  to 
the  inclusion  of  conductive  current  in  the  sourceless  frequency-domain 
Maxwell's  equations,  by  the  substitution  e>e  +  Tm*  In  the  language  of 
Faraday  lines  of  force,  this  means  that  every  E-fleld  line  also  becomes 
a  current-density  line.  Iftus,  both  capacitance  and  conductance  contain 
the  same  geometrical  factor,  X,  with  dimensions  of  length 

C  -  Xe 


G  =  Xo 


The  process  of  network  synthesis  is  considerably  simpler  for  a 
network  containing  only  two  element  types  (L-C,  R-C,  or  R-L)21  than  for 
a  network  containing  three  element  types  (L-R-C).  In  our  case,  although 
magnetic,  electric,  and  dissipative  mechanisms  are  all  present,  we  have, 
in  fact  (neglecting  radiation  resistance  for  the  moment),  only  two 
element  types: 

*  W.  Cauer,  Die  Verwirklichung  von  Wechselstroirwiderstanden 
vorgeschriebener  Frequenzabhangigkeit ,  Archiv  f.  Elektrotechnik,  17_ 
(1927),  355. 
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1.  Inductors — L  (Z  =  iaiL)  with  impedance  Z^*  iwl  and 
2*  "Repacitors"-- ( B)  with  admittance  Yg  =  =  iwc  —  G 

B  .  . 

=  iwc  1  -  — ^  . 

ea» 

(The  authors  have  coined  the  term  "repacitor"  to  denote  a  hybrid  circuit 
element  consisting  of  a  capacitor  and  resistor  connected  in  parallel. 
This  new  term  is  introduced  in  preference  to,  say,  "lossy  capacitor"  in 
order  to  emphasize  the  universality  of  the  medium -dependent  factor 
.  All  capacitors  become  conductive  according  to  the  same  ratio, 
as  long  as  the  material  parameters  are  spatially  homogeneous.)  Clearly, 
any  L-B  circuit  would  reduce,  in  the  case  of  zero  conductivity,  to  an  L- 
C  circuit,  and  in  the  case  of  high  conductivity,  to  an  L-R  circuit.  As 
long  as  the  medium  in  which  an  antenna  is  embedded  is  (1)  source  free 
and  (2)  time-invariant,  the  above  considerations  apply.  This  follows 
from  e  +  e  +  ■jjjj'  in  the  Maxwell  equations,  and  hence  should  correctly 
model  even  phenomena  arising  from  finite  skin  depth  in  the  high  o 
case.  Nor  should  the  inclusion  of  constant  resistors  (to  model 
radiation  resistance)  affect  the  validity  of  the  capaci tor-to-repacitor 
substitution.  Thus,  one  is  tempted  to  propose  the  following  algorithm 
for  generating  the  desired  circuit  model: 

(a)  Generate  (by  whatever  means)  the  h(oi)  and  Z(u>)  functions 
for  the  a*0  case. 

(b)  Synthesize  (by  whatever  means)  circuits  which  correctly 
model  the  h(u>)  and  Z(u>)  found  in  step  (a). 

(c)  Replace  all  capacitors  in  (b)  by  repacitors,  according  to 

the  formula  C  ♦  c(  1  -  . 

(d)  If  o  is  time-varying,  so  too  will  all  repacitors  become 
time-varying.  (At  this  point,  frequency-domain  concepts  lose  their 
validity. ) 

In  practice,  steps  (a)  and  (b)  are  very  vulnerable  to  error. 
It  is  important  to  realize  that  network  synthesis  is  not  a  unique 
process.  Thus,  a  circuit  which  correctly  models  the  o  =  0  case  need  not 
necessarily  model  the  more  general  a  *  0  case.  For  instance,  small 
errors  in  certain  capacitors,  which  might  not  be  apparent  in  the  non- 
dissipative  medium,  might  lead  to  large  errors  when  the  capacitors 
become  conductive.  Thus,  steps  (a)  and  (b)  should  be  repeated  for  a 
range  of  constant  conductivities  between  zero  and  the  peak  conductivi¬ 
ties  of  interest. 

Once  it  is  established  that  the  circuit  models  the  antenna 
accurately  for  the  time -invariant  case,  the  resistors  can  be  allowed  to 
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vary  in  time.  This  approach  has  yielded  very  good  agreement  both  with 
experimental  data  and  finite-difference  calculations  which  numerically 
solve  the  Maxwell  equations  in  the  neighborhood  of  the  antenna. 

In  effect,  allowing  a  to  vary  in  time  implies  that  there  are 
two  methods  of  driving  the  antenna: 

(1)  Directly,  through  the  incident  electromagnetic  field  and 

(2)  Parametrically,  through  o(t). 

The  second  method  is  analogous  to  the  behavior  of  a  parametric 
amplifier,  in  which  the  signal  modulates  a  capacitance.  Thus,  even  in 
an  electrostatic  field,  the  antenna  can  be  driven  by  o(t).  Physically, 
this  corresponds  to  a  conduction  current  source. 

Figures  24  to  27  show  the  effect  of  time-varying  air  conduc¬ 
tivity  on  antenna  coupling.  These  figures  show  measured  short-circuit 
current.  Figures  24,  25,  and  26  differ  in  the  relative  timing  of  AURORA 
and  the  pulser  firing  and  in  the  magnitude  to  o(t).  Figure  27  shows  a 
measurement  with  a  shorter  antenna. 
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Figure  24.  Comparison  of  measured  short-circuit  current  monopole 
antenna  response  (2.42  m)  with  and  without  ambient  time-varying  ionized 
air  conductivity  (time  in  s  x  1 0—  8 ) . 
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Figure  25.  Comparison  of  measured  short-circuit  current  monopole 


antenna  response  (2.42  m)  with  and  without  ambient  time-varying  ionized 
air  conductivity  (time  in  s  x  10”9). 


6.3  Theoretical  Interpretation 


The  experimental  results  have  been  interpreted  with  two  inde¬ 
pendent  calculational  approaches.  First,  the  short-circuit  current  was 
calculated  with  a  monopole  code  that  numerically  calculates  the  current 
response  from  the  incident  E-field  by  using  Maxwell's  equations  in 
f inite-di f ference  form  with  the  appropriate  boundary  conditions.  The 
code,  very  similar  to  that  developed  by  Merewether , 22  is  a  two- 
dimensional  code  that  takes  as  input  the  E-field  that  would  arise  in  the 
absence  of  the  antenna.  (This  formulation  has  been  widely  used  by  the 
EMP  community,  but  has  several  significant  disadvantages  for  systems 
applications.  First,  the  method  is  appropriate  for  the  solution  of 
simple  geometries,  but  for  realistic  Army  antennas,  it  can  be  quite 
intractable.  Similarly,  the  results  of  the  antenna  response  obtained  by 
this  method  are  difficult  to  interface  with  the  frequently  used  circuit 
analysis  codes — to  determine  the  system  response  at  the  component 
level.  This  method  has  obtained  respectable  results  for  simple 
geometries  in  underground  tests  and  therefore  is  used  here  as  a  basis 
for  judging  the  validity  of  the  equivalent  circuit  approach.) 
Comparisons  between  calculational  and  experimental  results  are  shown  in 
figure  28  (without  conductivity)  and  figure  29  (with  conductivity). 

Next,  calculations  were  made  using  the  LPN  equivalent  circuit 
obtained  by  the  algorithm  presented  in  section  6.2.  The  LPN  free-space 
antenna  LPN  equivalent  circuit  models  the  electrical  energy  and  magnetic 
energy  storage  and  the  dissipative  energy  loss  of  the  real  antenna. 
Theoretical  justification  of  this  approach  and  a  discussion  of  its 
limitations  are  presented  more  thoroughly  in  a  forthcoming  report. 
Assumptions  implicit  in  the  method  are  that  conductivity  is  independent 
of  field  intensity  (a  good  assumption  for  relative  humidity  of  about  40 
percent  and  fields  of  less  than  about  50  kv/m )  and  spatially  homog¬ 
eneous  . 


t 


2-2David  E.  Merewether,  Transient  Currents  Induced  on  a  Metallic  Body 
of  Revolution  by  an  Electromagnetic  Pulse,  IEEE  Trans.  EMC,  EMC-13 ,  2 
(May  1971). 
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Finally,  to  complete  the  cycle,  the  LPN  calculations  were 
compared  to  the  two-dimensional  finite  difference  procedure,  Compari- 
sions  were  made  over  a  wide  range  of  assumptions  about  conductivity 
waveforms.  Sensitivity  of  the  method  to  both  peak  conductivity  (always 
assuming  linearity,  i.e.,  =  0)  and  high-frequency  content  of  the 

conductivity  waveform  were  investigated.  Hie  comparisons  are  shown  in 
figure  30  without  conductivity  and  figure  31  with  'conductivi ty.  There 
was  some  doubt  about  the  applicability  of  the  method  beyond  the 
"adiabatic  regime,"  where  a(t)  varies  slowly  in  comparison  with  the 
driving  fields.  Figure  32  compares  Merewether’s  finite-difference 


Figure  30.  Comparison  between  LPN  and  finite  difference  code  response 
calculations.  No  time-varying  air  conductivity  present. 


Figure  31.  Comparison  between  LPN  and  finite  difference  code  response 
calculations.  Time-varying  air  conductivity  present. 


calculational  procedure  with  the  lumped -parameter  equivalent  circuit 
method  for  E(t)  and  o(t)  of  the  same  shape — a  distinctly  nonadiabatic 
condition.  Normalization  of  o(t)  was  chosen  so  that  the  peak  conduc¬ 
tivity  was  4  x  10-3  mho/m.  Agreement  is  remarkably  good. 

Thus,  a  three-way  agreement  has  been  established  between 
(1)  experiment,  (2)  a  finite-difference  method,  and  (3)  a  lumped- 
parameter  equivalent  circuit  method.  Experimental  results  have  also 
been  obtained  for  more  general  antennas — capacitively  and  inductively 
loaded  linear  antennas,  loop  antennas,  and  helical  antennas.  Hie 
authors  feel  that  the  lumped-parameter  method  can  be  successfully 
generalized  to  accommodate  these  configurations. 


7.  LIMITATIONS 


7. 1  Experimental  Apparatus 

Many  limitations  of  the  experimental  approach  have  been 
mentioned  in  the  previous  discussions  but,  for  the  sake  of  clarity,  they 
are  further  highlighted  here. 

First,  AURORA  produces  a  conductivity  time  history  whose  rise 
rate  is  slower  than  desired. 

The  parallel -plate  transmission  line  has  several  drawbacks  for 
this  type  of  application.  First,  the  volume  available  for  testing  (the 
region  between  the  plates)  is  limited,  since  the  line  must  be  placed  in 
the  AURORA  test  chamber.  The  plate  separation  cannot  be  too  great 
because  of  the  effect  of  the  test  chamber’s  conducting  walls. 
Furthermore,  in  order  to  successfully  launch  an  electromagnetic  wave 
down  the  line,  it  must  be  fed  from  a  point  through  a  tapered  transition 
section  before  it  reaches  the  parallel  section.  This  tapered  region 
cannot  be  used  for  testing,  yet  it  consumes  a  rather  large  section  of 
the  limited  available  space.  In  addition,  the  parallel-plate  line  must 
be  terminated  in  such  a  way  as  to  prevent  a  reflected  wave  from 
returning  to  the  paralleled  test  section.  A  conventional  way  of  doing 
this  is  by  tapering  the  line  and  installing  a  resistor  which  has  the 
value  of  the  characteristic  impedance  of  the  parallel-plate  line.  For 
our  experiment,  in  order  to  have  the  parallel  section  as  long  as 
possible,  the  termination  taper  section  was  omitted  and  the  transmission 
line  was  terminated  abruptly  in  an  "anechoic  curtain"  of  parallel  ana 
series  resistors.  Even  though  the  total  terminating  resistance 
approximated  that  of  the  line,  the  abrupt  geometry  change  resulted  in 
some  reflected  signal  (approximately  5  percent).  Although  not 
considered  significant,  its  presence  must  be  recognized. 

As  previously  discussed,  the  parallel-plate  line  itself  is 
affected  by  the  presence  of  the  conductivity,  which  tends  to  "short  out" 
the  field  in  the  line.  A  related  effect  that  is  present  has  been  called 
the  "inductive  kick."  This  is  a  sudden  rise  in  the  electric  field  after 
the  conductivity  is  removed.  It  is  present  because  of  the  inherent 
inductance  of  the  line.  This  "shorting  out"  and  "inductive  kick" 
behavior  is  a  manifestation  of  the  experimental  setup  and  not  really 
related  to  the  threat  environment.  However,  as  described  earlier,  the 
authors  are  confident  that  these  effects  can  be  greatly  reduced. 

Objects  placed  in  the  line  can  also  perturb  its  behavior.  For 
instance,  the  long  antenna  stretching  nearly  from  plate  to  plate 


distorts  the  current  in  the  hot  plate,  or  equivalently  induces  an  image 
antenna.23  We  have  established,  through  analysis  with  the  monopole 
code,  that  this  effect  is  negligible  for  the  cases  we  considered. 

Another  effect  which  is  present  and  must  be  addressed  is  the 
boundary  layer.  Since  electrons  are  not  easily  emitted  from  a  metal 
surface,  yet  are  highly  mobile  in  an  ionized  gas,  negatively  charged 
metallic  surfaces  exposed  to  a  radiation  environment  such  as  AURORA 
become  surrounded  by  an  electron  depletion  layer.  This  affects  the 
electromagnetic  coupling.  In  the  case  of  these  experiments,  the  effect 
can  be  ignored  because  of  properties  of  air  and  the  geometry.  However, 
in  the  event  that  experiments  dealing  with  smaller  geometries  such  as 
electronic  circuit  elements  or  some  other  medium  are  attempted,  the 
effect  may  be  pronounced. 

Another  thing  that  must  be  highlighted  concerns  the  nature  of 
the  experiments  that  have  been  conducted.  The  experiments  and  test 
apparatus  were  designed  to  consider  only  the  response  of  gamma-thin 
objects  (objects  which  are  thin  with  respect  to  a  gamma  path  length  and, 
therefore,  do  not  experience  any  direct  charge  buildup  due  to  the  radi¬ 
ation).  This,  then,  allows  one  to  make  the  approximation  that  the  en¬ 
vironment  and  the  coupling  are  separable;  that  is,  one  can  assume  that 
the  field  which  would  exist  without  the  system  (or  antenna)  present  is 
the  incident  field  which  is  then  coupled  to  the  system.  To  treat  the 
gamma  thick  problem24  (which  might  be  more  realistic  when  considering 
all  Army  systems)  would  be  considerably  more  difficult. 

Another  effect  which  is  present  in  the  actual  threat  situation 
is  that  of  neutrons.  The  neutrons  undergo  inelastic  scattering  as  well 
as  capture  and  result  in  further  late-time  sources  of  gamma  radiation. 
This  effect  has  been  totally  ignored.  This  assumption  is  probably  not 
too  severe  because  it  is  both  of  lower  magnitude  and  later  in  time  and 
probably  has  little  effect  on  the  short  antennas  considered  here — even 
though  it  must  be  considered  a  limitation  of  the  general  application  of 
this  scheme. 

When  viewed  in  the  context  of  experimental  verification  of 
analytical  models,  the  experimental  configuration  can  be  considered 
successful.  Hiis  is  not  to  say  that  the  LAEMP  environment  has  been 

23Clayborne  D,  Taylor,  George  A.  Steigerwald ,  On  the  Pulse  Excitation 
of  a  Cylinder  in  a  Parallel  Plate  Wave  Guide,  Mississippi  State 
University,  Sandia  Laboratories  Sensor  and  Simulation  Note  99,  Vol ,  7 
(March  1970), 

24G.  Merkel,  D,  J.  Spohn,  C.  L,  Longmire,  and  W.  F.  Crevier,  An 
Equivalent-Circuit  Analysis  of  the  HDL  Concentric  Cylinder  Experiments 
in  AURORA ,  IEEE  Trans  Nucl  Sci,  NS-2 4 ,  6  (December  1977), 
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simulated  in  the  sense  that  in  order  to  demonstrate  system  survivability 
one  need  only  expose  the  system  to  its  environment  and  then  demonstrate 
that  the  system  remains  functional.  This  was  not  the  goal  of  the 
experiment  and  should  not  be  interpreted  as  a  result. 

7.2  Equivalent  Circuit  Model 

There  are  several  well-known  limitations  of  equivalent-circuit 
models  (such  as  frequency  limitations)  which,  of  course,  apply  to  the 
formulation  presented  here.  In  addition  to  these  conventional  limita¬ 
tions,  this  model  demands  that  the  conductivity  be  spatially  homogeneous 
(even  though  it  could  be  generalized  to  nonhomogeneous  situations,  as 
illustrated  by  the  transmission-line  model  presented  in  section  4.1). 

In  a  related  issue,  the  coupling  to  antennas  is  naturally  a 
function  of  the  driving  field.  Uncertainties  in  the  environment  (par¬ 
ticularly  the  late-time  air  chemistry)  will  cast  a  similar  doubt  on  the 
antenna  coupling. 


8.  CONCLUSIONS  AND  EXPECTED  FOLLOW-ON  WORK 

The  authors  feel  that  considerable  progress  has  been  made  in  the 
understanding  of  source-region  EMP  coupling.  This  progress  has  taken 
the  form  of  advances  both  in  experimental  techniques  and  analytical  and 
numerical  modelling.  Certain  simplifying  assumptions  were  essential  to 
the  development  of  these  techniques.  As  is  often  required  in  research, 
many  simultaneous  effects  must  be  isolated  from  one  another  in  order  to 
gain  understanding. 

In  the  work  reported  here,  attention  was  concentrated  on  the  effect 
of  time-varying  air  conductivity  on  EMP  coupling.  Tactical  environ¬ 
mental  components  which  were  temporarily  neglected  in  both  experimental 
design  and  analysis  were  the  following. 

(1)  All  effects  of  local  Compton  currents  (system  charging,  space 
charging,  spatial  inhomogeneities  in  o(t>). 

(2)  All  neutron  effects. 

(3)  Boundary  layers  formed  at  metallic  surfaces  due  to  electron 
depletion  in  high  surface  fields. 

(4)  Avalanching  in  high  surface  fields. 

(5)  Nonlinearity  (E-field  dependence)  of  o(t). 


•  •  '  •  9 
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The  transmission-line  simulation  produced  excellent  results,  despite 
some  undesirable  effects,  i.e.,  line-AURORA  interaction  and  line-antenna 
interaction.  These  loading  effects  can  and  will  be  minimized  in  future 
experiments,  probably  using  techniques  described  above  (sect.  4). 

The  measurements  of  source-region  antenna  coupling  were  very 
instructive,  and  have  been  nicely  modelled  by  the  equivalent-circuit 
method  described  in  this  paper.  In  obtaining  experimental  data,  the 
authors  turned  to  AURORA  as  a  source  of  large-volume  ionizing  radiation, 
and  were  constrained  as  a  result  of  a  certain  conductivity  waveform 
whose  rise-rate  was  not  ideal  (too  slow).  However,  developments  in 
pulsed  power — at  HDL  (per  Alexander  Stewart)*  and  elsewhere — lead  us  to 
believe  that  better  sources  will  become  available.  In  any  case,  compar¬ 
isons  between  substantially  different  theoretical  methods  imply  that  our 
equivalent-circuit  treatment  can  be  generalized  to  the  case  of  faster- 
varying  conductivity. 

The  most  dramatic  effect  of  time-varying  conductivity  is  that 
conductive  current  (oE)  dominates  displacement  current  (c  in  driving 
an  antenna  with  a  low-impedance  load.  Whereas  the  displacement  current 
driver  tends  to  excite  resonances,  the  conduction  current  driver  tends 
to  be  low  frequency  and  monopolar. 

The  master-slave  transmission-line  approach  must  be  developed  for 
investigating  electromagnetic  coupling  in  time-varying  air  con¬ 
ductivity.  The  use  of  a  highly  ionizable  gas,  such  as  N2,  may  make  the 
technique  suitable  for  source-region  strategic  EMP  research. 

The  authors  also  expect  to  continue  to  study  the  behavior  of  ionized 
gases  at  various  pressures,  in  an  attempt  to  relate  macroscopic  (linear 
and  nonlinear  conductivity)  and  microscopic  (plasma  frequency,  collision 
frequency,  attachment  and  recombination  rates,  avalanching,  etc.)  param¬ 
eters.  Appropriate  experimental  approaches  may  include  capacitive 
probes,  "pie-pans,"  and  microwave  propagation  techniques.2^  For 
example,  nonlinear  conductivity  could  be  studied  by  looking  for 
frequency-doubling  in  microwave  signals. 

The  equivalent-circuit  antenna  analysis  method  will  be  applied  to 
more  general  antenna  types,  such  as  loops  and  log-periodic  antennal. 
The  technique  applies  whenever  o(t)  is  spatially  homogeneous. 
Techniques  for  spatially  inhomogeneous  environments  must  be  developed, 
if  gamma-thick  objects  in  the  source  region  are  to  be  treated. 

The  authors  are  developing  a  similar  series  of  code  verification 
experiments  for  communication  cables  exposed  to  the  threat  criteria 

25C.  E .  Baum,  Some  Limitations  on  Microwave  Air-Conductivity 
Measurements ,  Electromagnetic  Pulse  Sensor  and  Simulation  Notes,  Vol .  I, 
Note  16  (June  1970). 

*A.  Stewart,  Harry  Diamond  Laboratories,  private  communication. 


environment.  This  effort  presents  more  of  a  challenge  since  the 
orientation  of  the  field,  presence  of  the  ground,  and  lengths  involved 
all  must  be  treated  differently  than  the  simple  antenna.  Nevertheless, 
initial  indications  are  that  a  sufficiently  good  simulation  can  be 
accomplished  to  at  least  verify  the  coupling  codes  for  the  simplest 
configurations. 
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